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I.  MODELING  OF  IONIZED  GAS  SYSTEMS 
Program  Leader  -  Dr.  A.  V.  Phelps 

Our  studies  of  the  modeling  of  ionized  gas  systems  include  a)  the 
theoretical  prediction  of  the  conditions  for  the  growth  of  a  constriction 
in  a  weakly  ionized  gas,  such  as  found  in  a  molecular  gas  laser;  b)  the 
development  of  theories  of  radiation  transport  and  hydrodynamics  applic¬ 
able  to  a  wider  variety  of  real  situations,  such  as  expanding  and  turbu¬ 
lent  plasmas  in  geometries  other  than  the  usual  planar  geometry;  c)  the 
testing  of  radiative  transport  theories  under  controlled  laboratory  con¬ 
ditions;  and  d)  the  experimental  determination  of  the  laws  of  radiative 
scattering  and  absorption  appropriate  to  radiative  transport  problems. 

The  overall  objective  of  this  program  is  to  supply  experimentally  tested 
understanding  aud  computational  techniques  required  for  the  analysis  and 
prediction  of  the  role  of  a)  the  generation  and  transport  of  radiation 
in  plasmas  and  b)  the  transport  and  reactions  of  charged  particles  in  weak 
to  fully  ionized  gases.  Because  of  its  concern  with  a  number  of  simul¬ 
taneous  atomic  processes  and  with  relatively  large  numbers  of  interacting 
atoms  or  molecules,  this  program  also  provides  an  important  link  between 
the  experimental  and  theoretical  determinations  of  cross  sections  and  rate 
coefficients  which  make  up  much  of  the  remainder  of  the  program  and  the 
technology  of  ionized  and  excited  gases  which  this  work  is  designed  to 


support. 
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S^bility  of  Discharges  In  Weakly  Ionized  (Dr>  E>  F>  Jaegar> 

Dr.  L.  Osrer  and  Dr.  A.  V.  Phelps) 

The  objective  of  this  theoretical  project  is  the  quantitative 
Investigation  of  factors  which  control  the  growth  of  instabilities  in 
weakly  ionised  gas  discharges,  such  as  used  in  high  power  molecular 
lasers.  The  present  phase  of  this  project  is  concerned  with  instabilities 

which  result  from  a  decrease  in  gas  density  and  the  accompanying  increase 

In  the  ratio  of  the  elertrlr  j  *.  ..  * 

ctrlc  field  to  the  gas  density  E/N  and  in  the 

rate  of  ionization. 

The  charged  and  neutral  particle  equations  describing  the  growth 
of  ionisation  and  the  flow  of  charged  and  neutral  particles  during  the 
transient  following  the  application  of  voltage  to  a  weakly  ionised 
gas  have  been  formulated  and  programmed  for  solution  on  a  hign  speed 
computer.  In  their  present  form,  the  program  solves  four  simultaneous 
equations  of  which  three  describe  the  energy,  momentum  and  number 
balances  for  the  neutrals  and  one  describes  the  number  balance  for  the 
electrons  and  positive  ions.  Infinite  cylindrical  geometry  is  assumed 
and  the  axial  electric  field  is  allowed  to  vary  as  expected  for  a  constant 
external  resistance  in  series  with  a  fixed  voltage  supply.  Ihe  electron 
transport  and  Ionisation  rate  coefficients  are  assumed  to  be  determined 
by  the  ratio  ot  the  instantaneous  axial  electrical  field  to  the  local 
gas  density  as  expected  for  high  pressure,  molecular  gas  lasers. 

Solutions  obtained  thus  far  have  been  for  a  gas  with  the  properties  of 
helium  as  given  in  a  previous  investigation  of  glow  discharge  stability 
by  Ecker  et  al.  so  as  to  make  possible  comparisons  with  the  earlier 

(196 7) '  W'  KCS11’  SpatSchek  and  Zeller,  Phys.  Fluids  10,  1037 
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results.  More  realistic  transport  and  rate  coefficients  for  helium  and 
molecular  gas  mixtures  are  currently  being  generated  under  the  program 
discussed  in  Section  II  of  this  report  and  will  be  applied  to  this 
problem  during  the  next  report  period. 

Results  obtained  using  this  computer  program  for  helium  at  1000  torr 
are  illustrated  in  Figs.  1-2.  These  curves  show  the  effect  of  a  change 
in  the  radial  distribution  of  the  initial  ionization  and  the  effect  of  a 
change  in  the  external  series  resistance  on  the  growth  of  the  electron 
density  at  the  axis  of  the  discharge.  Figure  2  shows  representative 
spatial  distributions  for  the  electron  and  gas  densities.  Figures  1  and  2 
show  results  for  an  initial  field  strength  of  1500  V/cm,  a  series  resistance 
of  1000  ohm  per  cm  of  discharge  length  and  an  initial  ionization  distri¬ 
buted  according  to  a  Bessel  function  with  zero  density  at  the  tube 

10  -3 

wall  and  a  density  of  10  cm  on  axis.  Figure  1  shows  that  the  electron 
density  rises  to  a  nearly  constant  value  in  about  20  nsec.  The  magnitude  of 
this  quasi-steady  state  electron  density  is  determined  largely  by  a  balance 
between  the  ability  of  the  circuit  to  supply  current  and  the  steady  state 
current  density  vs  E/N  relation  characteristic  of  the  gas  under  study.  This 
steady  state  current  density  or  input  power  density  per  molecule  vs 
E/N  relation  can  be  calculated  for  any  weakly  ionized  gas  mixture  in 
which  the  electron  loss  by  diffusion  and  electron  production  by  cumu¬ 
lative  ionization  are  small  provided  the  appropriate  electron-ion  recombin¬ 
ation  coefficient  is  known.  The  electron  density  distribution  curves 
ne/ne(r=0)  of  Fig.  2  show  that  during  this  quasi-steady  state  period  the 
electron  density  becomes  relatively  independent  of  position.  A13o, 
the  n -utral  atom  censity  curves  (not  plotted  for  these  early  times)  show 


Axial  Electron  Density,  n  (r=0)  and  Atom  Density,  N(r=0) 


Spatial  variation  of  the  electron  and 
neutral  atom  density  at  various  times 
during  the  transients  of  Fig.  1. 
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relatively  little  change  from  their  initial  constant  value. 

After  the  quasi-steady  state  current  (about  0.5  A)  flows  for  about 
500  ysec  the  energy  input  (about  200  J/JL)  is  sufficient  to  raise  the  gas 
temperature  significantly.  The  resultant  outward  flow  of  gas  raises 
the  E/N  and  ionization  rate  near  the  axis  and  lowers  the  E/N  a:  d 
ionization  rate  near  the  wall.  The  discharge  thus  constricts  toward  the 
axis  as  shown  in  Fig.  2.  For  the  discharge  conditions  used  in  this 
example  the  discharge  continues  to  constrict  until  the  current  channel 
becomes  too  narrow  for  accurate  calculations.  The  curves  shown  for 
700  ysec  are  only  roughly  correct  because  of  the  use  of  a  large  artificial 
viscosity  term  inserted  into  the  equations  in  order  to  reduce  the  computer 
time  for  the  longer  runs. 

A  few  calculations  have  been  made  for  cases  in  which  the  initial 
ionization  is  restricted  to  a  channel  of  much  smaller  diameter  than  the 
distance  to  the  tube  wall  as  is  expected  for  TEA  lasers  with  pin  type 
cathodes.  The  general  behavior  of  the  discharge  is  similar  to  that  shown 
in  Fig.  2,  i.e.,  an  initial  expansion  of  the  discharge  to  about  twice 
the  initial  diameter  is  followed  by  a  collapse  to  a  narrow  channel. 

However,  when  all  other  parameters  are  unchanged  the  time  scale  for 
the  rise  in  gas  temperature  is  about  an  order  of  magnitude  shorter  when 
the  diameter  initial  ionization  is  reduced  by  an  order  of  magnitude. 

Although  the  total  energy  input  decreases  with  radius  there  is  some  tendency 
for  the  energy  input  per  unit  volume  prior  to  a  given  temperature  rise 
to  increase  with  decreasing  radius  of  initial  ionization  and  with  decreasing 
external  resistance. 

As  reported  in  the  Quarterly  Management  Report  for  October  1 
through  December  31,  1972  these  theoretical  studies  of  discharge  stability 
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have  led  us  to  propose  that  the  stability  of  a  pulsed  discharge  in 
C02  can  be  significantly  improved  by  the  addition  of  a  small  amount  of 
a  detaching  gas  such  as  H2  or  CO,  i.e.,  the  addition  of  a  gas  which 
detaches  the  electron  from  0  ions  formed  in  the  dissociative  attachment 
process  in  C02.  The  high  degree  of  sensitivity  of  a  niff use  discharge  in 
an  attaching  gas  to  the  growth  of  instabilities  is  the  result  of  the  very 
rapid  change  in  the  net  ionization  frequency  with  E/N  which  occurs  at  the 
E/N  appropriate  to  low  and  moderate  current  densities  in  a  recombination 
dominated  discharge  (see  above).  We  have  tested  this  hypothesis  by 
experiments  carried  out  under  another  contract.  Thus,  we  found  discharges 
in  mixtures  of  He,  C02  and  H2  in  the  ratios  of  90:9:1  to  be  much  more 
reproducible  than  were  discharges  in  90%  He  and  10%  C02>  These  measure¬ 
ments  were  carried  out  at  atmospheric  pressure  with  axisymmetric  Rogowski 
shaped  electrodes  spaced  at  about  2  cm  and  with  an  applied  voltage  of 
about  11  tcV/cm. 

All  of  the  preceding  numerical  results  should  be  considered  prelimi- 
nary  pending  further  checks  of  the  computer  code  and  more  extensive  varia¬ 
tion  of  the  discharge  parameters. 


— a-— er^n%  Transport  of  Resonance  Radiation  (Dr.  A.  V.  Phelps, 

Dr.  J.  E.  Jenkins  and  Dr.  R.  L.  Smick) 


The  measurements  of  the  collision  broadening  of  potassium  resonance 
lines  by  potassium  and  of  the  specular  reflection  at  the  potassium-glass 
interface  have  been  completed  and  were  presented  at  the  Conference  on 
Spectral  Lines.  Soon  afterwards  Dr.  Jenkins  returned  to  England  and  Dr. 
Smick  began  assembling  apparatus  for  the  measurement  of  the  spectral 
distribution  of  scattered  radiation  using  laser  excitation.  The  long- 
awaited  cw  dye  laser  was  delivered  by  the  manufacturer.  As  expected, 
the  line  width  of  the  laser  output  was  significantly  greater  than 
desired  for  our  experiments.  Accordingly,  some  effort  has  been 
devoted  to  the  development  of  a  fixed  etalon  for  use  inside  the  laser 
cavity  and  to  the  assembly  of  a  rapidly  swept,  tunable  etalon  for  monitor¬ 
ing  the  laser  operation.  At  present  we  are  able  to  obtain  reasonably 
stable  laser  operation  with  effective  line  width  of  about  0.01  ym.  This 
value  is  comparable  with  the  resolution  of  our  spectrometer  and  should  be 
satisfactory  for  our  present  needs.  A  new  oven  for  the  Na  vanor  cell 
has  been  built  with  windows  for  the  incident  and  reflected  light 
beams  and  for  observation  of  the  diffuse  scattering.  The  data  recording 
equipment  has  been  revised  to  make  it  compatible  with  an  available 
multichannel  scaler.  It  is  expected  that  this  system  will  be  operational 
during  the  next  report  period. 
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Radiation  Hydrodynamics  (Dr.  D.  C.  Hummer,  Dr.  J.  I.  Castor  and  Mr.  P. 

Kunasz) 

Una,  formation  in  static  spherical  tMMn.  A  vary  ganaral  and  effi- 
ciant  node  for  solving  lina  transfar  problems  in  „  extended  static  spharioal 
cloud  has  been  conpletad  and  extensively  tasted.  This  coda  utilises  the  non¬ 
linear  iteration  technique  (the  so-called  variable  Eddington  factor  method) 
used  previously  for  continuum  problems.1'2  generalised  to  the  case  of  non¬ 
gray  opacity  with  coupling  between  frequencies.  An  account  of  this  work  is 
being  prepared  for  publication. 

To  illustrate  the  effects  of  geometry  upon  the  radiation  field,  we 
present  some  results  for  a  rather  artificial  model  that  allows  comparison 
with  the  standard  plane-parallel  case.  We  consider  a  spherical  shell  with 
inner  radius  R.  and  outer  radius  R  surrounding  a  hollow  core.  The  mean 
optical  thickness  of  the  snail  is  T  and  the  electron  temperature  and  density 
are  constant.  The  opacity  is  regarded  as  a  general  function  of  radius  and 
frequency;  for  the  present  purpose  we  choose  k(v,r)  •  K*(v)r'2  ,  where  *(v) 
is  the  Doppler  profile  normalized  to  unity.  Figure  3  shows  how  the  radial 
behavior  of  the  source  function  (which  is  proportioned  to  yy  is  modified 
as  R/R.  increases  from  unity,  the  plane -parallel  limit,  in  an  atmosphere  of 
fixed  0£tical  thickness  T,  here  taken  to  be  1000.  Two  cases  e  5  (y/C  + 

A21  10  and  c  -  10  are  considered.  In  the  former  case,  in  which  the 

product  cT  »  i  ,  most  of  the  photons  created  are  destroyed  in  the  gas 

(energy  returned  to  electron  gas),  so  that  at  the  inner  parts  of  the  cloud 

1  mi!"'  !n'i  B’  EyblCk1,  M°nthly  N°tices  R°y'  Ast™-  Soc.  152,  1 
9-(1971)!Slnelii  °°d  G'  H™”er'  Monthly  Notices  Roy.  Astron.  Soc.  153, 


H 


the  excitation  is  essentially  independent  of  geometry.  In  the  outer  part 
of  the  cloud,  the  excitation  decreases  steadily  as  the  size  of  the  cloud 
increases.  Fc  e  =  10  ,  eT  <<  1,  so  that  most  of  the  photons  escape. 

In  this  case,  the  level  of  excitation  throughout  the  cloud  depends  on 
geometry  quite  strongly. 

In  the  special  case  of  a  homogeneous  sphere  (constant  density) ,  the 
integral  for  the  source  functions  can  be  transformed  into  a  form  that  is 
amenable  to  solution  by  the  kernel  approximation  method."1  This  calculation 
is  very  easy  and  straightforward  to  program,  and  gives  results  that  should 
be  accurate  to  better  than  one  percent.  While  these  results  were  first 
obtained  in  order  to  check  the  more  general  code,  they  are  of  some  interest 
in  themselves,  for  a  comparison  of  the  source  functions  with  those  for  plane- 
parallel  slabs  whose  optical  thickness  is  equal  to  the  optical  diameter  of 
the  sphere  show  the  following  result: 

SsPh(T)  "  CSpp<T> 

when  C  <  1  and  is  practically  independent  of  T,  i.e.  depends  only  one,  T 
and  the  frequency-dependence  of  the  opacity.  The  smallest  value  of  C  encoun¬ 
tered  in  a  very  exhaustive  survey  of  the  parameter  space  for  both  Doppler  and 
Lorentz  profiles  was  ~  0.44.  A  similar  result  was  obtained  by  Avery, 

House  and  Skumanich  for  finite  cylinders  from  Monte  Carlo  calculations. 

Some  light  on  this  result  is  obtained  from  the  consideration  of  the  pro— 

bability  that  an  excited  atom  will  emit  a  photon  that  escapes  the  cloud  with  no 

3 

E.  H.  Avrett  and  D.  G.  Hummer,  Monthly  Notices  Roy.  Astron.Soc.  130,  294 
C 1965) « 

4 

L.  W.  Avery,  L.  L.  House  and  A.  Skumanich,  J.  Quant.  Spectrosc.  Radiat. 

Transfer  519  (1969). 
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further  scattering.  We  have  shown,  fcr  any  distribution  of  density  with 
radius,  that  the  probability  P^Ct)  is  related  to  the  corresponding  quantity 
for  the  plane-parallel  case  for  P^Ct),  ^y  t*ie  relation 


P  .  (t)  =  P  (t)Q(r) 
sph  pp 

where  Q(r)  is  an  integral  containing  the  radial-dependence  of  the  density 
and  satisfies  the  following  inequalities: 


1  5  Q(r)  ^  2  ,  Doppler  profile 


and 

1  ^  Q(r)  ^  3/2  ,  Lorentz  profile. 

Thus  the  escape  probability  depends  primarily  on  the  shortest  path  to  the 
surface  and  only  very  weakly  on  longer  paths.  It  appears  that  this  circum¬ 
stance  might  provide  the  basis  for  a  systematic  calculational  procedure  for 
handling  radiative  transfer  problems  in  complicated  geometries. 

For  the  results  shown  in  Fig.  3  in  which  the  level  of  excitation  does 
depend  strongly  on  the  size  of  the  sphere,  some  other  factor  must  be  operative. 
This  seems  to  be  a  dilution  effect,  since  the  power  input  goes  as  R  (for 
k(r)  ~  r-2),  while  the  available  volume  goes  as  R3.  This  point  is  now  being 

examined  in  detail. 

line  formation  in  expanding  spherical  geometries-  While  the  variable- 
Eddington  factor  procedure  can  be  generalized  to  handle  radiative  transfer 
spherical  atmospheres  with  velocity  fields,  a  number  of  problems  have  to  be 
overcome.  In  order  to  obtain  accurate  solutions  that  can  be  used  to  check 
the  results  from  the  variable-Eddington  factor  procedure,  we  have  devised  a 
less  general  non-iterative  technique  that  can  treat  velocities  up  to  a  few 
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times  larger  than  the  mean  thermal  velocity.  This  method  depends  a 
generalization  of  Rybicki's'  rearrangement  of  the  standard  Feautrier6  proce¬ 
dure  for  solving  the  discrete  form  of  the  radiative  transfer  equations.  Al¬ 
though  the  matrix  of  the  system  is  very  large,  all  but  one  of  the  submatrices 
are  tri-diagonal  or  have  some  other  useful  special  form,  so  that  the  amount  of 
computing  required  is  surprisingly  small.  This  procedure  has  been  coded  and 
is  now  nearly  debugged.  In  the  limit  of  zero  velocity,  it  should  yield  a 

good  check  on  the  Eddington-f actor  results  for  the  static  case,  since  it  is 
a  completely  independent  procedure. 


.  G*  B*  Rybicki’  J*  Quant*  Spectrosc.  Radiat.  Transfer  11,  589  (1971). 
P.  Feautrier,  C.  v.  hebd.  Seanc.  Acad.  Sci.  Paris  258,  3189  (1964). 
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Plasma  Statistics  (Dr.  W.  E.  Brittin  and  Dr.  w.  R.  Chappell) 

Work  was  done  on  the  problem  of  formulation  of  fundamental  theories 
designed  to  treat  interaotlng  composite  particles,  such  as  electrons,  protons, 
atoms,  ions,  etc.  Although  no  completely  satisfactory  theory  was  discovered 
rough  this  work,  it  was  felt  that  this  work  provided  part  of  the  basis  for 
subsequent  work  which  appears  to  be  satisfactory.  Basically,  the  general 
idea  was  an  attempt  to  treat  various  composite  particles  as  though  they 
were  elementary  particles  with  internal  structuie.  In  essence  what  was 
sought  was  a  generalised  quantum  mechanical  transformation  which  would 
allow  one  to  pass  from  the  standard  representation  which  describes  the  elemen¬ 
tary  particles  comprising  the  system  to  a  new  basis  described  in  tents  of 
composite  particle  quantities. 
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II.  TRANSPORT  COEFFICIENTS  AND  REACTION  RATES 
Program  Leader  -  Dr.  S.  J.  Smith 

The  prediction  of  the  performance  characteristics  of  ionized  gas 
devices,  such  as  the  efficient,,  power  output  and  stability  of  high  power 
lasers,  requires  that  one  be  able  to  calculate  accurate  values  of  elec¬ 
tron  transport  and  reaction  rate  coefficients  for  electrons  and  ions.  In 
order  to  meet  some  of  these  requirements  we  are  engaged  in:  a)  a- determin¬ 
ation  of  sets  of  electron  collision  cross  sections  consistent  with  measured 
electron  transport  and  ionization  coefficients  and  with  recent  electron 
beam  measurements  of  excitation  cross  sections  for  the  upper  excited  states 
of  gases  found  in  molecular  lasers;  b)  measurement  of  cross  sections  for 
molecular  dissociation;  c)  evaluation  and  compilation  of  low  energy  electron 
photon,  and  ion  collision  cross  section  data;  and  d)  determination  of  nega¬ 
tive  ion  characteristics,  such  as  single  and  multiple  photodetachment  cross 
sections,  negative  ion  stability  and  low  electron  energy  scattering  para¬ 


meters  . 
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Electron  Transport  and  Ionization  Coefficients  (Dr.  G.  E.  Chamberlain, 

Dr.  L.  J.  Kieffer  and  Dr.  A.  v.  Phelps) 

For  the  past  six  months  we  have  continued  in  our  effort  to  summarize 
and  analyze  the  techniques  for  measuring  electron  transport  and  ionization 
coefficients.  We  are  compiling  a  complete  set  of  these  types  of  data  for 
N2  in  order  to: 

(a)  compare  the  data  taken  by  differing  techniques 

(b)  look  for  systematic  errors 

(c)  be  able  to  test  against  the  data  our  criteria  for 
data  selection  and  estimates  of  systematic  errors 
as  these  criteria  and  estimates  develop. 

We  have  concentrated  on  measurements  in  Nj  since  we  believe  they  exemplify 
most  of  the  problems  that  occur  in  the  measurement  of  transport  and  ioniza¬ 
tion  coefficients  and  because  of  the  importance  of  N2  in  molecular  gas  dis¬ 
charge  lasers.  We  have  put  together,  in  first  draft  form,  several  compo¬ 
nents  of  our  review. 

(A)  A  summary  of  the  principal  techniques  (drift  tube  with  electrical 
shutters  and  without  shutters)  for  measuring  drift  velocity  and  a  dir  cuss ion 
of  the  corrections  to  the  primary  data  that  are  necessary  to  account  for 
effects  of  diffusion,  ionization,  absorbing  boundaries,  and  high  gas  densities. 

(B)  A  summary  of  the  Townsend-Huxley  technique  for  measuring  the 
characteristic  energy  (diffusion/mobility)  by  observing  the  radial  spread 

of  a  swarm.  The  high  sensitivity  of  the  method  of  geometric  alignment,  space 
charge,  and  attaching  gas  impurities  are  noted.  To  this  is  to  be  added  a 
discussion  of  systematic  effects  in  the  Townsend-Huxley  method  due  to  the 
assumed  boundary  conditions.  An  additional  section  is  to  be  written  on 
diffusion/mobility  values  obtained  from  drift  gap  measurements,  although  the 
technique  involved  is  already  covered  in  the  discussion  of  drift  velocities 
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obtainea  by  the  drift  gap  technique. 

(C)  We  have  begun  a  discussion  and  wade  some  estimates  of  the  effects 
of  space  charge  with  and  without  ionization  on  the  measured  values  of  trans¬ 
port  and  ionization  coefficients.  Drift  velocities  are  at  least  susceptible 
to  space  charge;  no  experimental  evidence  has  been  found  for  a  change  in 
measured  drift  velocity  due  to  high  electron  densities.  Ionization  and  dif¬ 
fusion  measurements  are  certainly  sensitive  to  space  charge  and  a  simple 
model  has  been  considered.  Additional  material  in  the  literature  exists  and 
will  be  added  to  the  space  charge  review.  Our  current  efforts  are  directed 
toward  completing  the  revxew  of  diffusion  measurements,  initiating  and 
completing  a  review  of  ionization  measurements,  and  augmenting  the  review 
and  modeling  of  space  charge  effects. 

Electron  transport  and  rate  coefficient  bibliography  scheme.  The  index¬ 
ing  scheme  for  maintaining  the  Information  Center's  bibliography  of  data  on 
electron  transport  and  rate  coefficients  is  defined  by  the  list  of  index 
terms  and  descriptors  shown  in  Table  1.  The  bibliography  will  be  sorted  on 
the  ii.dex  terms,  the  descriptors  will  be  attached  to  each  citation  of  the 
paper  as  supplementary  information  on  the  type  of  paper  and  type  of  measurement. 

The  data  categories  are  meant  to  be  general  categories  indicating  the 
type  of  coefficient.  For  example,  measurements  of  (diffusion  coefficient)/ 
(mobility),  or  characteristic  energy,  are  to  be  found  under  Diffusion.  The 
Collision  Frequencies  category  has  been  added  to  cover  the  numerous  types  of 
elastic  and  inelastic  collision  frequencies  that  occur  in  the  literature. 

Our  coverage  is  expanded  to  include  energy  relaxation  measurements  which  will 
be  fou.  i  under  Excitation,  or  Collision  Frequencies  as  appropriate  to  the 
specific  paper. 

We  nave  dropped  the  Conductivity  category.  Papers  with  measurements  of 
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conductivity  ho.  which  drift  veiocity  cr  collision  fluency  values  cm  „e 
obtained  will  be  indexed  in  the  appropriate  category. 

The  Analysis  descriptor  has  been  added  to  identify  papers  in  which 
electron  collision  cross  sections  have  been  determined  by  making  comparisons 
with  electron  transport  and  rate  coefficient  data. 


Table  1 

Index  Terms  and  Descriptors  for  Bibliography  of  Data  on 
Electron  Swarm-Gas  Interactions 

INDEX  TERMS 


(1)  Data  Categories 

Drift  Velocity 

Diffusion 

Excitation 

Ionization 

Attachment 

Detachment 

Recombination 

Collision  Frequencies 

Distribution  Function 

(2)  Species 

Atoms 

Molecules 

Mixtures 

(3)  Citation  Number 


descriptors 

Experimental 

Theoretical 

Analysis 

Drift 

Afterglow 

Discharge 

Breakdown 
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JILA  Information  Analysis  Center  (Dr.  L.  J.  Kieffer) 

In  addition  to  our  continuing  to  keep  current  bibliographies  on  low 
energy  electron  collision  cross  sections  and  photoabsorption  and  ionization 
cross  sections  and  compiling  the  data  for  these  processes,  we  have  initiated 
a  new  project  on  electron  transport  and  rate  coefficients.  A  current  bibli¬ 
ography  and  data  compilation  will  be  maintained.  This  is  a  follow-up  to  a 
project  which  was  carried  on  for  many  years  by  Professor  J.  Dutton  and  which 
is  being  terminated  by  a  critical  data  compilation  mentioned  in  previous 
reports.  We  have  devised  a  new  indexing  scheme  (see  Table  1)  and  expanded 
the  coverage.  This  project  is  being  carried  out  by  Dr.  G.  E.  Chamberlain 
of  the  Information  Center  Staff. 

We  continue  to  help  Dr.  A.  V.  Phelps  by  preparing  best  cross  section 
sets  for  discharge  laser  modeling.  Some  revisions  in  the  currently  used 
nitrogen  cross  section  have  been  made  but  more  work  will  be  done. 

Work  has  begun  in  preparing  a  set  of  atomic  hydrogen  cross  sections 
for  use  in  modeling  electrical  discharges  in  I^. 

Molecular  Dissociation  Processes  (Dr.  L.  J.  Kieffer) 

No  experimental  work  was  done  under  this  project  during  this  reporting 
period.  Data  taken  last  July  (see  previous  semiannual  report)  were  analyzed 
Two  manuscripts  are  being  prepared  for  publication. 


groups  we 


_Electron  Transport  in  Mercurv  Vapor  (Dr.  A.  V.  Phelps) 

As  the  result  of  numerous  inquiries  from  laser  development 
have  assembled  and  reviewed  most  of  the  available  data  on  electron 
transport  and  ionization  coefficients  in  mercury  vapor.  We  have  also 
reviewed  many  of  the  available  references  dealing  with  the  spectra 
emitted  by  the  first  few  excited  states  of  mercury  and  with  the  collision 
processes  leading  to  the  production  of  these  states. 

The  available  electron  transport  coefficient  data  for  mercury 
consist  primarily  of  measurements  of  electron  mobility  and  electron 
temperature"  as  determined  from  measurements  in  the  positive 
columns  of  low  pressure  discharges.  There  is  only  one  measurement 
of  electron  drift  velocities  using  the  drift  tube  technique1  and 
this  author’s  results  show  a  large  scatter  from  one  graph  to  another. 

Two  groups  have  measured  Townsend  ionization  coefficients,  but  there 
appears  to  be  no  measurement  of  the  ratio  of  the  diffusion  to  mobility 
coefficients.  It  is  strongly  recommended  that  any  attempts  to  determine 
a  set  of  electron  collision  cross  sections  consistent  with  the  positive 
column  data  be  carried  out  with  a  computer  code  capable  of  including  the 
effects  of  electron-electron  scattering  on  the  electron  energy  distribution. 
Unfortunately,  our  code  is  for  weakly  ionized  gases  only  and  does  not 
include  the  effects  of  electron-electron  scattering. 

At  present  there  seems  to  be  no  completely  consistent  model  of 
radiation  sources  and  collision  kinetic  which  will  describe  the  available 
data  on  fluorescence  in  mercury  vapor.  It  is  therefore  recommended 
that  appropriate  spectroscopic  and  kinetic  experiments  be  undertaken. 

1  E.  W.  McCutchen,  Phys.  Rev.  112.  1848  (1958). 
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Three  experiments  which  seem  particularly  valuable  are  a)  the  determination 

O 

of  the  molecular  corstants  of  the  states  of  Hg2  formed  from  the  Hg  P 
3 

and  Pq  levels  of  Hg  using  laser  induced  fluorescence  of  higher  Hg^ 

2 

states  as  proposed  by  Hessel  and  Smith,  b)  the  measurement  of  far  wing 

fluorescence  following  optical  excitation  of  the  P  state  and  c)  the 

measurement  of  far  wing  absorption  at  high  mercury  temperatures,  i ,  e. .  , 

3 

a  modern  version  of  the  Kuhn  and  ^reudenberg  experiment.  With  the 
results  of  this  combination  of  experiments  one  should  be  able  to 
construct  a  set  of  potential  energy  curves  and  transition  probabilities 
appropriate  to  the  excited  states  of  mercury  of  laser  interest. 


M.  M.  Hessel  and  E.  W,  Smith,  proposal  to  Air  Force  Weapons  Laboratory. 
H.  Kuhn  and  K.  Freudenberg,  Z.  Physik  76,  38  (1932). 
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Negative  Ion  Kinetics  (Dr.  W.  C.  Lleeberger.  Dr.  J.  Hall,  Dr.  H.  Hotop,  and 


Mr.  T  A.  Patterson) 


During  the  present  contract  period  we  have  completed  our  studies  of 
metallic  native  ions  of  gold,  platinum,  copper,  and  silver.  TVo  manu¬ 
scripts  have  been  prepared1’2  which  describe  these  results;  these  papers 
will  be  published  in  the  Journal  of  Chemical  Physics.  Further  studies  of 
metallic  negative  ions  using  our  new  sputtering  ion  source  will  be  pursued 
when  the  laser  photoelectron  spectroscopy  apparatus  is  available. 

The  studies  of  the  electron  affinities  of  the  noble  metals  suggest  that 
we  should  be  able  to  obtain  relatively  accurate  estimates  of  the  elements  in 
the  three  long  series:  (K-Cu,  Rb+Ag,  Cs-Au)  by  means  of  horizontal  iso- 
electronic  extrapolations.  In  particular,  Charkin  and  Dyatkina3  and  Zollweg4 
have  shown  that  the  energy  difference  between  the  two  configurations  dV  and 

k 

J  s  (k=0,...  ,10)  increases  nearly  linearly  with  k.  Since  it  is  well  estab¬ 
lished  that  the  ground  state  configuration  of  negative  ions  in  the  three  long 
periods  is  d  s2,  one  can  estimate  the  electron  affinities  of  these  atoms  if 
those  for  K,  Kb,  Cs  on  the  one  side  and  those  for  Cu,  Ag,  and  Au  on  the  other 

side  are  known,  because  thereby  the  slope  of  the  straight  line  representing 
k  2  k 

the  d  s  +d  s  energy  differences  between  the  negative  ion  and  the  dks  neutral 
atom  state  is  fixed.  The  results  of  such  an  extrapolation  are  presented 
in  Table  2,  and  probably  represent  the  best  present  estimates  of  the  elec¬ 
tron  affinities  of  the  elements  in  the  three  long  series. 

For  this  extrapolation  we  had  relied  upon  the  alkali  electron  affinities 

H.  Hotop,  R.  A.  Bennett  and  W.  C.  Lineberger,  J.  Chem.  Phys.  (March  1,  1973). 

2 

H.  Hotop  and  W.  C.  Lineberger,  J.  Chem.  Phys.  (March  1,  1973). 

3 

0.  P.  Charkin  and  M.  E.  Dyatkina,  J.  Struct.  Chem.  6_,  397  (1965). 

4 

R.  F.  Zollweg,  J.  Chem.  Phys.  50,  4251  (1969). 


wents  in  the  three  long  aeries  [eV] 


calculated  by  Weiss.5  While  we  would  like  to  measure  the  electron  affini¬ 
ties  of  the  alkali  metals  with  our  tunable  laser  apparatus  it  is  not  current¬ 
ly  possible  for  us  to  make  2  ym  tunable  light,  which  would  be  required  to 
see  the  photcdetachment  thresholds.  However,  recent  calculations  by 
Norcross  and  Moores  have  shown  that,  at  the  thresholds  for  the  first  excited 
state  channels  in  the  alkalis,  there  should  be  a  pronounced  cusp-like  feature. 
The  location  of  this  cusp  would  thus  provide  an  accurate  measurement  of  the 
electron  affinities  of  the  alkalis.  Further  these  cusps  will  appear  in  an 
energy  region  accessible  to  our  tunable  laser  apparatus.  Consequently,  we 
have  constructed  a  double  charge  transfer  alkali  negative  ion  source.  Our 
prelimir.ary  measurements  on  sodium  show  that  this  cusp  is  in  fact  present 
and  imply  that  EA(Na)  =  0.542  +  0.010  eV.  Similar  measurements  on  other 
alkalis  are  currently  in  progress. 

During  the  present  reporting  period,  we  have  been  very  fortunate  to 
have  several  JILA  visitors  working  on  the  megawatt  pulsed  nitrogen  laser. 
Professor  V.  Ehlers  has  completed  the  testing  and  performance  evaluation  of 
the  laser  and  a  manuscript  describing  the  laser  will  soon  be  submitted  to  the 
Review  of  Scientific  Instruments.  Dr.  R.  Keller  has  successf' 1 ly  used  the 
nitrogen  laser  to  pump  dyes  which  laze  in  -he  8,000-9,000  A  region.  By  use  of 
several  dyes  as  excitation  transfer  agents,  he  has  been  able  to  achieve  fairly 
efficient  conversion  (greater  than  5%)  from  3371  A  to  8500  A.  This  infrared 
laser  will  be  used  to  determine  the  electron  affinity  of  atomic  oxygen. 

While  waiting  for  machine  shop  modifications  to  the  new  charge  transfer 
source  system,  some  existing  data  have  been  analyzed  to  yield  electron 
affinity  information  for  the  ions  ,  NH  ,  OH  ,  SC^  ,  and  S2  •  The  resul- 

5  A.  W.  Weiss,  Phys.  Rev.  122,  1826  (1961);  166,  70  (1968). 

^  E.  Clementi,  Phys.  Rev.  136 ,  A980  (1964). 
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tant  affinity  values  cannot  lay  claim  to  such  high  precision  and  absolute 
accuracy  as  our  careful  measurements  and  analysis  of  NO-  and  02".  However, 
the  errors  associated  with  the  experimental  data  are  only  £  10  meV,  the 
primary  uncertainties  being  in  the  interpretation  of  rotational  effects. 

All  of  the  data  have  been  taken  relative  to  0_  and  have  been  expressed  in 
absolute  units  using  the  0  result  of  Hotop,  Patterson  and  Lineberger.  Thus 
we  report  EA(NH~)  =  (0.779  +  0.037:>eV,  EA(OH)  =  1.829^,°^  eV, 

EA(NH)  =  (0.38  +  0.03)eV,  EA(S02)  =  (1.097  +  0.036)eV,  and  EA(S2)  = 

(1.663  +  0.040)eV. 

Additionally  the  angular  distributions  were  measured  for  OH-  and  NH2” 
to  yield  the  anisotropy  parameter  8  at  4880  A  as  6(NH  ”)  =  0.027  ±  0.012 
and  8 (OH  )  =  -0.993  +  0.040.  Our  data  show  no  0H(l,0)OH~  peak  in  the  photo¬ 
detachment  spectrum  larger  than  0.0006  of  the  strong  peak  OH(0,0)OH-.  Thus, 
using  arguments  given  by  Branscomb7  we  can  calculate  an  upper  limit  on  the 

change  in  intemuclear  separation  upon  photodetachment.  We  find  Ar  £  3  x 

e 

-4  « 

10  A  formally,  although  the  theoretical  model  probably  loses  validity 
before  this  level  is  reached.  This  work  is  in  the  final  stages  of  preparation 
for  publication. 

During  the  next  report  period  we  hope  to  survey  a  large  number  of 
interesting  atomic  negative  ion  photodetachment  spectra,  using  the  new 
"Coulatron"  universal  positive  ion  source  and  a  charge  transfer  system  which 
should  soon  be  operational. 


L.  M.  Branscomb,  Phys.  Rev.  166 ,  11  (1966). 


7 
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III.  ENERGY  LOSS  AND  RADIATION  PRODUCTION  BY  ELECTRON  COLLISIONS 
Program  Leader  -  Dr.  E.  C.  Beaty 

Interaction  of  Electrons  with  Ions  (Dr.  G.  H.  Dunn  and  Dr.  F.  L.  Walls) 

The  measurement  of  the  energy  dependence  of  electron-molecular  ion 
recombination  cross  sections  is  being  carried  out  with  an  ion  trap  tech¬ 
nique  illustrated  in  Fig.  4.  During  this  report  period  a  new  trap  fabricated 
from  molybdenum  to  high  tolerances  was  installed  and  tested.  Figure  5  shows 
an  ior.  peak  of  mass  18  (NH^) ,  and  demonstrates  achievable  mass  resolution 
and  signal  to  noise  with  the  new  trap.  Indications  are  that  ion  lifetimes 
and  signal-to-noise  ratios  are  now  adequate  to  pursue  recombination  measure¬ 
ments  as  planned.  Since  September  1,  1972  this  work  has  been  supported  by 
another  agency. 


NOISE  POWER 


Fig.  5.  Noise  power  spectrum  showing  trapped  ion  resonance  at  f  =  o)z/2tt 
for  mass  number  18  (NH4+) .  Distinction  from  H20+  has  been  made 
by  cyclotron  resonance. 


-29- 


Theoretical  Atomic  Physics  (Dr.  D.  G.  Hummer  and  Dr.  D.  W.  Norcross) 

Electron  excitation  cross  sections  for  atmospheric  ions.  in  this 

period,  a  considerable  amount  of  effort  was  expended  in  improving  the  atomic 

structure  code  to  allow  more  accurate  atomic  models  to  be  obtained.  The 

wave  functions  for  C  III,  N  IV  and  0  V  are  significantly  better  than  those 

obtaireu  in  our  early  calculations  and  seem  to  be  comparable  to  those  obtained 

by  other  very  sophisticated  procedures.  The  calculated  oscillator  strengths 

agree  quite  well  with  the  results  of  beam-foil  spectroscopy  (generally  to 
within  30%). 


The  Distorted  Wave  Code  is  also  being  modified  to  provide  better  accu¬ 
racy,  greater' convenience  in  use  and  substantially  increased  efficiency,  in 
preparation  for  the  larger  calculations  that  we  are  planning  for  the  immediate 
future.  We  expect  to  begin  production  within  the  next  two  months. 
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Atomic  and  Radiation  Interaction  Theory  (Dr.  Sydney  Celt  ran) 

Free-free  absorption  coe ff icienrs ■  We  have  completed  the  work  on 
the  calculation  of  atomic  free-free  absorption  coefficients  and  It  has  been 
accepted  for  publication.1  The  coefficients  have  been  calculated  for  the 
electron-neutral  atom  systems  involving  He.  C.  »,  0,  He,  Ar,  Kr,  and  Xe. 
Some  angular  distributions  are  presented  and  thermal  averages  have  been 
evaluated  in  the  ranges  of  X  -  .5  to  20  pm  and  T  *  500  to  20000°K.  The 

computer  codes  remain  available  for  application'  to  any  neutral  atom  of 
interest. 

Angular  distributions  for  Ionization.  The  "Coulomb-projected  Bom" 
approximation  is  being  applied  to  the  high  energy  electron  impact  Ionization 
of  the  hydrogen  and  helium  atoms.  This  method  has  proved  useful'in  describ¬ 


ing  the  large-angle  differential  cross  sections  for  proton-hydrogen  charge 
transfer2  and  electron  Impact  excitation  of  hydrogen  and  helium,3  where  the 
conventional  Born  approximation  breaks  down  completely.  The  mathematical 
problems  of  evaluating  the  ionization  amplitude  are  much  more  severe  than  in 
the  case  of  excitation,  but  progress  is  being  made. 

Interaction  of  atoms  with  intense  radiation  fields.  The  usual  form  of 


perturbation  theory  which  is  applied  to  the  calculation  of  absorption  coeffi¬ 
cients  breaks  down  at  radiation  intensities  of  the  order  of  1014  w/cm2, 
where  the  interaction  energy  between  electron  and  field  is  comparable  to  the 
binding  energy  between  electron  and  nucleus.  Since  laser  intensities  of 
order  1019  w/cm2  are  now  routinely  available,  it  is  urgent  that  the  absorption 
nschanism  ofjuch  intense  radiation  be  thoroughly  understood.  We  are  inves- 

1  S.  Geltman ,  J.  Quant.  Spectrosc.  Radiat.  Transfer  (in  press). 

2 

S.  Geltman,  J.  Phys.  B  4_,  1288  (1971). 

3 

S.  Geltman  and  M.  B.  Hidalgo,  J.  Phys.  B  4,  1299  (1971);  5,  617  (1972). 
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tisating  alternate  theoretical  procedures  in  which  the  electron-nuclear 
interaction  is  regarded  as  the  perturbation,  while  the  electron-field  inter. 


action  is  taken  to  be  part  of  the  unperturbed  problem.  The  method  is  being 


applied  to  the  N-photon 
hydrogenic  transitions, 
cal  approaches.^ 


absorption  in  bound-bound,  bound-free  and  free-free 
and  the  results  will  be  compared  with  other  theoreti- 


H.  Reiss,  Phys.  Rev.  A  1. 
(1972). 


803  (1970);  G.  J.  Pert,  J. 


Phys.  A  5,  1221 


IV.  GENERATION  AND  ABSORPTION  OF  RADIATION 
Program  leader  -  Dr.  A.  C.  Gallagher 


Gas  Line  Broadening  and  Intermolecular  Potentials  (Dr.  A.  C. 

Gallagher  and  Dr.  C.  G.  Carrington) 

We  have  measured  the  pressure  dependence  of  the  Rb-Xe  molecular- 
continuum  radiation  to  determine  the  rate  of  formation  of  bound  Rb*-Xe 

molecules : 

*  * 

Rb  +  2Xe  -*•  RbXe  +  Xe. 

This  step  following  electron  excitation  of  the  atom  is  an  important 
part  of  the  association  laser  dynamics.  Combined  with  the  equivalent  measure¬ 
ments  of  the  formation  rates  for  Na*  and  Li*  that  are  planned,  this  should 
allow  us  to  understand  how  this  process  depends  on  well  depth,  vibrational 
spacing,  and  other  parameters  of  possible  association  laser  molecules. 

Quantum  mechanical  calculations  of  Frank-Condon  factors  are  being 
carried  out  to  determine  the  validity  of  the  classical  model  we  have  proposed 
for  calculating  the  optical  properties  (e.g.,  gain)  of  association  laser 
molecules.  Preliminary  results  for  the  bound-free  and  bound-bound  transi¬ 
tions  are  shown  in  Fig.  6.  In  the  spectral  region  of  primary  interest  for 
association  lasers  this  indicates  about  10%  accuracy  of  the  classical 
approximation.  In  some  other  regions  such  as  the  far  red  edge  of  the  spec¬ 
trum  there  are  major  differences. 
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Classical  vs  quantum-mechanical  emission  spectrum  for  radiatively  dissociating  mole 
cules.  The  potentials  are  modeled  after  the  Cs-Ar  system;  the  temperature  is  half 
the  well  depth.  The  bound  AB*  curves  represent  the  portion  of  the  total  radiation 
due  to  bound  excited  molecules.  The  relation  between  the  spectrum  and  the  AB*  poten 
tial  is  shown  in  the  insert. 
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Excitation  of  Metal  Resonance  Lines  (Dr.  A.  C.  Gallagher  and  Mr.  D.  Leep) 

We  have  made  preliminary  measurements  of  the  excitation  cross  section 
of  the  Li  resonance  line.  The  Li-rare  gas  molecules  are  expected  to  have 
the  best  interaction  potentials  of  the  various  alkalis  for  use  as  an  associa¬ 
tion  laser  molecule.  Such  a  laser  would  utilize  electron  excitation  of  the 
lithium  resonance  line,  with  efficiency  dependent  on  this  cross  section. 
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